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Abstract 

We address the problem of classification of integrable differential-difference equations in 
2+1 dimensions with one/two discrete variables. Om" approach is based on the method 
of hydrodynamic reductions and its generalisation to dispersive equations as proposed in 
[lOl [H] . We obtain a number of classification results of scalar integrable equations including 
that of the intermediate long wave and Toda type. 
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1 Introduction 



This paper is aimed at the classification of scalar 2 + 1 dimensional integrable equations of the 
general form 

ut = F{u,w), (1) 

where u{x, y, t) is a scalar field, w{x, y, t) is the nonlocal variable, and F is a differential/difference 
operator in the independent variables x and y (the explicit form of w and F will be specified in 
what follows). Publications |im lllj provide a novel perturbative approach to the classification 
of integrable equations of the form ([1]) which have non-degenerate dispersionless limit. This 
approach is based on the requirement that all hydrodynamic reductions [S] of the dispersionless 
limit can be deformed into reductions of the full dispersive equation (a brief description of this 
method is included in Sect. 2). For the simplest nonlocality w = D~^DyU (or, equivalently, 
Wx = Uy), the paper [llj gives a complete list of integrable third order PDEs of the form 

Ut = (pux + ipUy + r]Wy + e(...) + e^(...), (2) 

where ip, and rj are functions of u and w, while the terms at e and are assumed to be 
homogeneous differential polynomials of the order two and three in the x— and y— derivatives 
of u and w. 

Theorem 1 ^ll] Up to invertible transformations, the examples below provide a complete list 
of integrable third order equations (0j with non- degenerate dispersionless limit: 

KP equation ut 
modified KP equation ut 

Gardner equation ut 
VN equation ut 

modified VN equation ut 

HD equation ut 

deformed HD equation ut 
Equation Ei ut 
Equation E2 ut 
here B = fSuD^ — Dy, j3, 5=const. 

Although most of these examples are well-known [17^ [3T] , the last three equations are apparently 
new. We refer to [TT] for technical details, Lax pairs, etc. Based on the same method, in [23] the 
above classification was extended to fifth order scalar equations, while the paper jl6j provides a 
complete list of integrable two-component systems of Davey-Stewartson type (all with the same 
nonlocality Wx = Uy). See also [H [191 [26] ^^'i references therein for the existing classification 
results of integrable systems in 2+1 D. 
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Remark. The search for integrable equations ^ with 'nested' non-locahties of the form Wx = 
Uy, Vx = f{u,w)y leads to commuting flows of equations from Theorem 1, and does not give 
essentially new examples. 

The aim of this work is to extend the above classification to various differential-difference 
non-localities w such as 

A^w = ^ ^ — Uy, Wx = AyU, A^u; = u. 

The first two cases are known as the intermediate long wave and the Toda type non-localities, 
respectively. We use the following standard notation for the e-shift operators and the for- 
ward/backward discrete derivatives, 

Txf{x,y) = f{x + e,y), f{x,y) = f{x - e,y), A+ = A" = i^l^, 

same for Ty, Ty^, A+, Ay. Note that all our non-localities reduce to Wx = Uy in the disper- 
sionless limit e — >• 0. 

Here is a brief summary of our classification results (see Sect. 3 for the corresponding Lax 
pairs, dispersionless limits, etc). In Sect. 3.1 we consider two classes of equations with the 
nonlocality of intermediate long wave type. Our results are summarised in Theorems [2] and [3] 
below: 



Theorem 2 The following examples constitute a complete list of integrable equations of the 
form 

Ut = ifUx + ipUy + TWx + iJWy + e(...) + e^(...), 

where w is the non-locality of intermediate long wave type, A^w = "^^^^Uy, the coefficients 
ip,ip,T,r] are functions ofu,w, and the terms at e and are assumed to be homogeneous differ- 
ential polynomials of the order two and three in the x— and y— derivatives of u and w: 

Ut = UUy -\- Wy, 

Ut = {w + ae^)uy + Wy, 

Ut = U^Uy + {Uw)y + —Uyyy, 

2 . ^ 

Ut=U Uy+ {Uw)y + -^\Uyy 

The first example appeared in |7] as a differential-difference analogue of the KP equation, see 
also |27J. It can be viewed as a 2 -|- 1 dimensional integrable version of the intermediate long 
wave equation [35]. The last two examples are differential-difference versions of the Veselov- 
Novikov and the modified Veselov-Novikov equations, respectively (the third example appeared 
previously in [25]). The second and the fourth equations seem to be new. 

Theorem 3 The following examples constitute a complete list of integrable equations of the 
form 

Ut = ipuy + rjwy + fA^g pA~q, 
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where w is the non-locality of intermediate long wave type, A^w = ^^y^Uy, and il^,ri, f, g,p,q 
are functions of u and w: 

Ut = UUy + Wy, 

Ut = {w + ae^)uy + Wy, 

Ut = wuy + Wy + — ^e^"*, 
Ut = wUy + Wy + e"(A+ + A;^)e". 

Here the first two examples are the same as in Theorem 2, while the third appeared in [IB]. The 
fourth equation is apparently new. 

In Sect. 3.2 we consider the case of integrable equations with the Toda non-locality. The 
main result is as follows: 



Theorem 4 The following examples constitute a complete list of integrable equations of the 
form 

Ut = ifUx + fA^g+pAyQ, 
where w is the non-locality of Toda type, Wx = AyU, and (p,f,g,p,q are functions of u and w: 

Ut = uAyW, 

Ut = (au + /3)A-e"', 

Ut = e'"^A+^ + ^A-{e'"V^), 

here a, /3 = const. 

The first example is the well-known Toda equation, while the second is equivalent to the Volterra 
(Toda) chain when a 7^ (a = 0), respectively. The third equation seems to be new. 

In Sect. 3.3 we discuss the case of integrable equations with two discrete variables. Our 
main result is 



Theorem 5 The following examples constitute a complete list of integrable equations of the 
form 

Ut = fA+g + hA-k + pA+q + rAyS, 

where w is the fully discrete non-locality of the form A^w = AyU, and f,g,h,k,p,q,r,s are 
functions of u and w: 

Ut = uAy{u - w), 

Ut = u{A^ -\- Ay)w, 

Ut = (ae-" + /3) A-e"-'", 

nt = (ae" + /3)(A+ + A-)e"', 

Ut = - ;5e2« (^e"'-"A+^a-/3e2« + Ay {e'^-^'y^a - /3e^^)'^ , 
here a, (3 = const. 
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In equivalent form, the last example is known as the 2 + 1 dimensional analogue of the modified 
Volterra lattice [30]. 

Although fully discrete 3D equations are outside the scope of this paper, they can be dealt 
with in a similar way, suggesting an alternative to the standard approach based on the multidi- 
mensional consistency [21 [3l [23l [29] . The only constraint required by our method is that of the 
non-degeneracy of the dispersionless limit (see Sect. 2). Numerous examples of this type are 
provided by various equivalent forms of the Hirota equation |15[ I22j governing KP/BKP/Toda 
hierarchies. These include the following difference equations (we prefer to write them using 
discrete derivatives A rather than shifts T: in this form the dispersionless limit becomes more 
easily seen): 

(AJ u - AJu)A+ A+t, + (A+u - AJ u)A+A+« + (AJu - A+i,)A+ A+i, = 0, 
Aj(lnfi::)+Aj(l.ffj)+A,+ (l„ftH)=0, 

a e^-^y"" + /3 e^-^t^ + 7 e^^^t^ = 0, 
a e^- " + f3 e^y^y"" + 7 e^^^*"" = 0, 

etc, here q,/3,7 are arbitrary constants. The corresponding dispersionless limits e — t- can be 
obtained by replacing discrete derivatives by partial derivatives. This gives 

{Ux - Uy)Uxy + {Ut - Ux)Uxt + {Uy - Ut)Uyt = 0, 




a e"-2' + /3 e"-' + 7 e^f* = 0, 



a e"-- + /3 e'^yy + 7 e"" = 0, 

respectively. In different context, these and similar dispersionless equations appeared in [ll[5l[6l 
113] . see also references therein. One can show that all of the above difference equations inherit 
hydrodynamic reductions of their dispersionless limits, at least to the order e^. This property 
can be used to classify integrable difference equations within particularly interesting subclasses. 
As a simple illustration let us consider equations of the form 

AtA-{u) + a; Ay if {u)) + AtA;{g{u)) = 0, 

where / and g are functions to be determined. The requirement that hydrodynamic reductions 
of its dispersionless limit, Uxx + f{u)yy + g{y)tt = 0, can be deformed into reductions of the 
difference equation up to the order e^, leads to the following constraints for / and g: 

f" + g" = 0, g"{l + f')-g'f" = 0, /"2(i + 2/')-/'(/' + l)r = 0. 

Without any loss of generality this gives f{u) = ln(l — e") — u, g{u) = — ln(l — e"), resulting in 
the difference equation 

AtA-{u) + A+A-(ln(l - e«) - u) - A+A,-(ln(l - e")) = 0, 



5 



which is yet another equivalent form of the Hirota equation, known as the 'gauge-invariant 
form' [32] , Its dispersionless Umit appeared recently in the classification of integrable equations 
possessing the 'central quadric ansatz' [T^. We hope to report further classification results 
elsewhere. 



2 Classification scheme 

Let us illustrate our approach using the example of the 2+1 dimensional Toda equation, 

Ut = uAyW, Wx = A+n. 
Expanding the right hand sides using Taylor's formula one obtains 

Ut e 

— = — -Wyy + —Wyyy + . . . , 

-If! 

Wx — Uy + ^Uyy + ^ Uyyy + . . . . 



(3) 



The corresponding dispersionless limit results upon setting e = 0: 

Ut = UWy, Wx = Uy. (4) 

This dispersionless system admits exact solutions of the form 

u = R, w = w{R), (5) 

where R{x, y, t) satisfies the pair of Hopf-type equations, 

Ry = fiRx, Rt = jJ^RRx- (6) 

Here is an arbitrary function, and w' = fi. Solutions of this type are known as one-phase 
solutions (or planar simple waves, or one-component hydrodynamic reductions). One can show 
that both relations ([6]) can be deformed into solutions of the full Toda equation in the form 

u = R, w = w{R) + ewiRx + e'^{w2Rxx + W3Rl) + 0{€^), (7) 

and 

Ry =iiRx + e^{aiRxxx + ct2RxRxx + a^Rl) + O(e^), 
Rt =^l^RRx + e^iPiRxxx + hRxRxx + PsRl) + O(e^), 

where Wi,ai,/3i are certain functions of R. We point out that, modulo the Miura group [H], 
the relation u = R can be assumed undeformed. Furthermore, one can show that odd order 
e-corrections in the equations ([8]) (but not d?])) must vanish identically. Substituting ([7]) into ([3]), 
using ([8]) and the compatibility condition Ryt = Rty, one can explicitly calculate all coefficients 
in d?! and ([8]) in terms of /u and its derivatives by collecting terms at different powers of e |10] . 
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This gives 



^ 2 

W3 ~ [R {^i'f (2// - V) + (11^' + V)) , 

/32 =^Rfi (r {fi'Y - 2Rfx') + V (3/x' + Rfi")) , 

k =^R (i? {n'f (2/i - + SRfi^n'n" + II {ii'f - 3i?V) + /^^ (V + R^J'^^^)) , 

etc. We point out that this calculation is an entirely algebraic procedure. Thus, all one-phase 
solutions of the dispersionless system are 'inherited' by the original dispersive equation, at least 
to the order O(e^): it is still an open problem to prove their inheritance to all orders in the 
deformation parameter e. It is important that this works for arbitrary n{R). The requirement of 
the inheritance of all hydrodynamic reductions of the dispersionless limit by the full dispersive 
equation is very restrictive (even to the order 0(e^)), and can be used as an efficient classification 
criterion in the search for integrable equations. In all examples considered so far, the existence 
of such deformations to the order was already sufficient for integrability (in many cases, even 
the order was enough), and implied the existence of conventional Lax pairs. 

As an illustration, let us consider the class of Toda-type equations of the form 

ut = f^yd, = A+n 

where f{u,w) and g{u,w) are two arbitrary functions. Imposing the requirement that all one- 
phase solutions of the corresponding dispersionless limit are inherited by the full dispersive 
equation, we obtain very strong constraints for / and g. Indeed at the order e we get 

9u ~ 0) fufw ~ Oj fw if9ww "I" 9wfw) ~ 0, 

so that, excluding the degenerate case /„ = 0, we arrive at gu = 0, /«, = 0. At the order one 
obtains two additional constraints: 

/» = 0, g"{wf-g'iw)g"iw) = 0. 

Modulo elementary changes of variables this leads to the cases f{u) = u, g{w) = w and f{u) = 
au + (3, g{w) = which correspond to the Toda and Volterra chains, respectively (see Section 
3.2). Note that all these constraints appear at the order e^, and are already sufficient for the 
integrability, implying the existence of Lax pairs. 
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Remark. All equations discussed in this paper possess dispersionless limits of the form 

Ut = ifUx + IpUy + r]Wy, Wx=Uy, (9) 

which will be assumed non- degenerate in the following sense: 

(i) The coefficient rj is nonzero: this is equivalent to the requirement that the corresponding 
dispersion relation defines an irreducible conic. 

(ii) The dispersionless limit ([9]) is not totally linearly degenerate. Recall that totally linearly 
degenerate systems are characterised by the relations 

rjw = 0, ilj^+r]u = 0, (Pw + ipu = 0, ipu = 0. 

Dispersive deformations of degenerate systems do not inherit hydrodynamic reductions, and 
require a different approach which is beyond the scope of this paper. We point out that most 
of the integrable examples of interest are non-degenerate in the above sense, or can be brought 
into a non-degenerate form. 



3 Differential-difference equations 

In this section we search for integrable examples within various classes of differential-difference 
PDEs generalising intermediate long wave and Toda type equations. We skip the details of 
calculations, which essentially follow the pattern outlined in Sect. 2. Classification results are 
presented modulo Galilean transformations, and transformations of the form u — ?■ au + /?, — t- 
aw + 7. 



3.1 Equations with the intermediate long wave non-locahty: A^w = 



2 



First we classify integrable equations of the form 

Ut = ifUx + ipuy + TWx + r]Wy + e(...) + e^(...), (10) 
where w is the non-locality of the intermediate long wave type, A+u) = ^y^Uy, or, equivalently. 



w 



f ^"^1 %- Here dots denote terms which are homogeneous polynomials of degree two and 
three in the x- and y-derivatives of u and w, whose coefficients are allowed to be functions of u 
and w. One can show that all e-terms, as well as all terms containing derivatives with respect 
to X, in particular ip and r, must vanish identically. 

Theorem 2 The following examples constitute a complete list of integrable equations of the form 
[W\) with the non-locality of intermediate long wave type: 

Ut = UUy+Wy, (11) 

Ut = {w + ae^)uy + Wy, (12) 

Ut = U^Uy+{UW)y + —Uyyy, (l3) 
Ut = U^Uy+{uw)y + —iUyy- --^\ . (l4) 
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Although the first two equations are of the first order, they should be viewed as dispersive: the 
dispersion is contained in the equation for nonlocality. The equation ()lip . which can be written 
in the form 

eTx + l 

Ut = UUy + -jr—^Uyy, 

first appeared in [7J as a differential-difference analogue of the KP equation, see also [27j (we 
point out that its dispersionless limit does not coincide with that of KP). It can also be viewed as 
a 2 + 1 dimensional integrable version of the intermediate long wave equation |35j . The equation 
()13p is a differential-difference version of the Veselov-Novikov equation discussed in |25) . The 
last example can be viewed as a differential-difference version of the modified Veselov-Novikov 
equation. To the best of our knowledge, equations (fT^ and (fH|) are new. 

Lax pairs, dispersionless limits and dispersionless Lax pairs for the equations from Theorem[2] 
are provided in the table below (note that equations and have coinciding dispersionless 
limits/dispersionless Lax pairs). Here and in what follows. Lax pairs were obtained by the 
quantisation of dispersionless Lax pairs as discussed in [M] . 



Eqn 


Lax pair 


Dispersionless 
limit 


Dispersionless 
Lax pair 


m 


Txip = etpy - uip 

2 

etpt = Y'^yy + (^ - f ■"s/)'^ 


Ut = UUy + Wy 
Wx = Uy 


e'^^ = Sy-u 
St = ^S^ + w 


m 


Txip = ee'^ipy - aip 

Ipt = ^Ipyy + {w - ^Uy)tl'y 


Ut = {w + ae'^)uy + Wy 

Wx = Uy 


^ e~^Sy - a 
St = ISy + WSy 


m 


e{Tx - l)iJy = -2u{Tx + 1)^ 
i't = ii'yyy + {w - ^Uy)iPy 


Ut = U^Uy + {Uw)y 
Wx = Uy 


eSa^+lSy - 2n 
St = 12 ^y ~^ wSy 


m 


e{Tx - l)Vy = f ^(T, - 1)V - 2n(r, + 1)^- 

= ^i'yyy + {W - ^Uy)lPy + ^{Wy - ^Uyy)tP 


Ut = U^Uy + {Uw)y 
Wx = Uy 


it^ , -, Sii 2li 

St = \2^y ujSy 



Remark. Equations (jlip and ()12p are related by a (rather non-trivial) gauge transformation. 
Let us begin with the dispersionless limit of ()12p . 

Ut = {W + ae^)Uy + Wy, W^ = Uy, 

with the corresponding Lax pair 

St = ^S^ + wSy, e^- = e-"5j, - a. 

Let h he a potential such that u = h^, w = hy. One can verify that the new variables u = 
w + ae", w = ht — S = S + h satisfy the dispersionless equation (fTT]l . 

Ut = UUy + Wy, Wx = Uy, 
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along with the corresponding Lax pair 



St = ^Sy + W, 6^=" = Sy-U, 



thus estabhshing the required hnk at the dispersionless level (it is sufficient to perform this 
calculation at the level of Lax pairs: the equations for u,w will be automatic). The dispersive 
version of this construction is as follows. We take the equation (jl2p . 

2 - 1 ' 



ut = {w + ae^)uy + Wy, w = --^ -Uy, 



with the corresponding Lax pair 

'^t = ^ipyy + {w - ^Uy)'il;y, T^ii = ee""V?; - oti). 

Let -ff be a potential such that u = -^^^i?, w = ^^2^ Hy. One can verify that the new variables 
u = Hy + oe", w = Ht - ^ + ^e^'A+Hy, -0 = e^/^tp satisfy the equation (fTTD . 



2 ^^y 

Ut = uuy + Wy, w = - ^ _i ^v^ 

with the corresponding Lax pair 

. £2 _ e . , , , 

Again, it is sufficient to perform this calculation at the level of Lax pairs. Due to the complexity 
of this transformation we prefer to keep both equations in the list of Theorem [2] as separate 
cases. 



Another interesting class of equations with the non-locality of intermediate long wave type 

is 

ut = ipuy + r]Wy + fA^g + pA~q, (15) 
where A^w = '^'^^^ Uy, and ip,!], f, g,p, q are functions of u and w. 

Theorem 3 The following examples constitute a complete list of integrahle equations of the form 
( f.?5|i with the non-locality of intermediate long wave type: 

Ut = UUy + Wy, 

ut = {w + ae^)uy + Wy, 

Ut = wny + Wy + ^^-^^e^"", (16) 
Ut = wUy + Wy + e''{A'^ + A'^)e''. (17) 

Here the first two equations are the same as in Theorem [2l the third example first appeared in 
|18| . while the fourth is apparently new. Lax pairs, dispersionless limits and dispersionless Lax 
pairs for equations from Theorem [3] are provided in the table below (note that equations ()16p 
and p7p have coinciding dispersionless limits): 
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Eqn 


Lax pair 


Disper sionless 
limit 


Disper sionless 
Lax pair 


m 




Ut = 2e'^^Ux + WUy + Wy 
Wx = Uy 


Sy = 2e" cosh Sx 
St = e2«smh25^+ 
2we'^ cosh Sx 


m 


we^{T^i) + T-^iP) + f e"[(A+ + A-)e«]V' 


Ut = 2e'^^Ux + WUy + Wy 
Wx = Uy 


Sy = 2e" cosh Sx 
St = e2"smh25^+ 
2we'^ cosh Sx 



3.2 Equations with the Toda non-locality: Wx = A+m 

In this section we classify integrable equations of the form 

ut = ^ux + fA+g + pAyq, (18) 

where the non-locahty w is defined as Wx = AyU, and 99, f,g,p, q are functions of u and w. 

Theorem 4 The following examples constitute a complete list of integrable equations of the form 
with the non-locality of Toda type: 

Ut = uAyW, (19) 
Ut = {au + ^)A-e'", (20) 
Ut = e'"V^A+V^+^A-(e"'V^), (21) 

here a, /3 = const. 

Equation (|19p is the 2+1 dimensional Toda equation, which can also be written in the form 
(lnu)xt = A+A~u, while equation (pO|) is equivalent to the Volterra chain when q / 0, or to 
the Toda chain when a = 0. Lax pairs, dispersionless limits and dispersionless Lax pairs for the 
equations from Theorem U] are provided in the table below: 
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Eqn 


Lax pair 


Dispersionless 
limit 


Dispersionless 
Lax pair 




CI in — 1 lln 

t-Lyipx — atp 

eipt = -Tytp + {T-^w)iIj 


ill — '^y 

Wx = Uy 


e Ox — 

St = -e^y + w 


& 


eTyipx = - iTyu)Tytp + {aTyU + /3)ip 
etPt = -e'^Tyi) + ae^V 


Ut = {au + /3)e'^Wy 

Wx = Uy 


e^ySx = -ue^y +au + fi 
St = -e^e^y + ae"" 


m 


eV^j = ie-T,V-i(T-ie-)r-V 


Ut = e^Uy + ue^Wy 

Wx = Uy 


e^y Sx = Sx - ue^y - u 
St = e^ sinh 5^^ 



Remark. One can show that there exist no non-degenerate integrable equations of the form 

Ut = rjWy + + pA~q, 

where the non-locahty w is defined as A+w = Uy, and r], f, g,p, q are functions of u and w. In- 
deed, the integrabihty requirement imphes the condition ij = 0, which corresponds to degenerate 
systems. 



3.3 Equations with the fully discrete non-locality: A^w = A+u 

In this last section we classify integrable equations of the form 

Ut = /A+5 + hA-k+ pA+q + rAyS, (22) 

where the non-locality w is defined as A^w = AyU, and the functions /, g, h, k,p, q, r, s depend 
on u and w. 

Theorem 5 The following examples constitute a complete list of integrable equations of the form 
\2S\) with the fully discrete non-locality: 



Ut = uAy{u — w), (23) 

Ut = n(A+ + A-)w;, (24) 

Ut = (ae-" + /3)A-e"-"', (25) 

Ut = (ae"+/3)(A+ + A;)e'", (26) 

Ut = - /3e2« (^e'^-^A+Va - /3e2« + A-(e"'-" y^a - /3e2«)) , (27) 



here a, (3 = const. 

In equivalent form, equation (j27|) is known as the 2 + 1 dimensional analogue of the modi- 
fied Volterra lattice [30]. Lax pairs, dispersionless limits and dispersionless Lax pairs for the 
equations from Theorem [5] are provided in the table below: 
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Eqn 


Lax pair 


Disper sionless 
limit 


Disper sionless 
Lax pair 




eipt = Tyip - wip 


Ut = u{Uy - Wy) 
Wx = Uy 


^Sx+Sy _ _^Sy _I^^^Sx 

St = e^y -w 


m 


T^Tylj) = Ty^l) - UIP 

eipt = Tyip + {Ty'^w)ip 


Ut = u{Uy + Wy) 
Wx = Uy 


^Sx+Sy ^ ^Sy _ ^ 

St = e^y +w 






ut = {a + I3e^)e~^{uy - Wy) 

Wx = Uy 


e-^-St = -e-'^St- 
ae-'"e-^- + ^e'"" 




Ty Ip — Q,e''+/3^a:V 1 ae^+p^ 
ni(T p'^^ih 


Ut = {ae"" + I3)e"{uy + Wy) 

WX = Uy 


^ ae"+/3 

e^-St = e-'^St- 

QpW _ w 
jje e Lxe 


m 


T^TyiP = ^{Tye-^)Ty'^+ 

= -ae'^Tyi, + l3{T-^e'")T-^i> 


Ui = a{e^--)y - P{e^+^)y 

Wx = Uy 


^Sx~\~Sy — — g ^Q^v 1 

St = -ae'"e^y + 
Pe'^e-^y 



Remark 1. The continuum limit of the modified Volterra lattice ()27p in x— direction, namely 
X — )• hx, u ^ hu and /i — )• 0, gives the Toda-type lattice (|2ip . Similarly, in the same limit 
equations ()23p and ()24p give the Toda equation ()19p . while the remaining two, ()25p and (|26p. 
lead to the equation ()20p with a = 0. 

Remark 2. We point out that there exist other types of integrable equations with the non- 
locality A+if = Ay li, which are not covered by Theorem [5l One of such examples is the first 
flow of the discrete modified Veselov-Novikov hierarchy constructed in [33] , 

Ut = Y^-(ry-iA+e2-)(A+e-2-), A+u; = A+n. 
This equation is not of the form (|22p . furthermore, its dispersionless limit is degenerate: 

Ut = 2Wx, Wx = Uy. 

4 Concluding remarks 

This paper outlines an approach to the classification of integrable differential-difference equations 
in 2+lD (with one/two discrete variables) based on the method of hydrodynamic reductions. 
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1. It would be challenging to extend our approach to the classification of fully discrete 3D 
equations. Some work in this direction, based on the concept of multidimensional consistency, 
has been done in [21 [29] , 

2. A novel approach to the integrability of multidimensional lattices of the Toda and Volterra 
type, based on the concept of characteristic Lie rings, was proposed recently in [T3]. Its relation 
to our method is yet to be properly investigated. 

3. Our calculations demonstrate that, even though nontrivial combinations of differential and 
difference operators (in the same independent variable) were allowed in the initial ansatz, non of 
such 'differential-delay' cases survived the integrability test, leading to either purely discrete, or 
purely differential situations. It remains to be seen whether there exist proper differential-delay 
integrable equations in 2+1 D. 
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